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ABSTRACT 

We explore molecular cloud properties and the physics of CO transition lines in z > 6 
Lyman-break galaxies and predict their CO fluxes using an analytic formalism built 
from global models of star formation in high-redshift galaxies that minimizes our 
reliance on local observations. Our model includes a new approach to calculating the 
molecular gas fraction that takes the total gas density, the star formation rate, and the 
star formation efficiency in clouds as the principle inputs. This method agrees with 
chemical equilibrium calculations of the molecular fraction based on local chemistry 
if galaxies at z > 6 have metallicities in the range of 3-30 percent of solar. Such low 
metallicitics in turn imply that much of the carbon in these systems exists in ionized 
form rather than as CO. Moreover, we find that CO transitions will typically be sub- 
thermally populated, especially those higher-order lines observable at high redshift 
with ALMA. Ultimately, we expect current facilities will not be able to observe the 
CO signal from reionization epoch galaxies except with great difficulty. We estimate 
that ~ 1000 hours of integration with JVLA will be required to detect the CO(l-O) 
transition in even the most UV-bright systems at z = 6. 

Key words: galaxies: high-redshift — galaxies: evolution — galaxies: ISM — ISM: 
molecules — radio lines: galaxies 



1 INTRODUCTION 



The 
has 



WFC3 camera 
discovered 



9 



aboard the Hubble Space Telescope 
large popula t ions of gala xies out to 



Bouwcns ct al.' '2006'; 'Bunke r et al] |2010|: 
McLure et al.l i2010,: Finkclstcin ct al.. .2010. : .Bouwens et al.l 
2011al lbl: lEllis et aLlfioTi T via the Lyman-break technique 



and opened investigation into galaxy formation during the 
first ~ 500 Myrs of cosmic time. Observations in the rest- 
frame optical and UV have been important in estimating 
star formation rates and stellar masses for these Lyman- 
break galaxies (LBGs) in addition to the census statistics of 
abundance and clustering. According to such measurements, 
typical galaxies during the epoch of reionization at z > 6 are 
much smaller and denser than th ose at low redshift (e.g., 
lOesch et al.ll2010l : IOno et al.ll2012h and would be considered 
dwarfs based on their stellar masses (e.g.. lLabbe et aLll2010l ). 
However, these observations tell us little about the thermal 
and chemical state of high-redshift interstellar media (ISM). 

On the other hand, high-resolution observations of 
molecular transition lines expected from the the Jansky Very 
Large Array (JVLA) and the Atacama Large Millimeter Ar- 
ray (ALMA), particularly those of CO and CII, will probe in- 
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ternal galactic physics and dynamics fe.g. lCarilli et al]|2008l : 
I Walter fc Carillill2008l '). Averaged over large scales, intensity 
maps of these lines can als o validate and enhance 21 cm 
studies of reioni z ation (e.g. iLidz et al] l201ll : ICarillil 1201 ll : 
iGong et al.ll201ll . l2012h . 

In this Paper, we probe the relationship between the 
conditions expected in the high-redshift ISM and the physics 
of CO transitions using t he model for high-redshift galac- 
tic disks we developed in iMufioz fc Furlanettol (|2012D . Ad- 
ditionally, we aim to make predictions for the CO signal from 
these sources to facilitate current and future observations. 

Recent attempts at predicting the CO emission sig- 
nal fr o m high-redshift galaxies (e.g., lObreschkow et ahl 
l2009bl : iLidz et al] l201ll ) have relied heavily on empirical 
relations calibrate d from local observations. For example, 
lObreschkow et all (l2009bl') used a semi-analytic model of 
galax y formation (|Croton et al.l l2006l : iDe Lucia fc BlaizotI 
l2007t ) to consider galaxy sizes and masses appropriate to 
high-redshift while including the effect of the cosmic mi- 
crowave background (CMB), whose temperature is nearly 
20 K by 2 = 6, as both a source of heating and as an observa- 
tional background. However, the galaxy formation prescrip- 
tions of their semi-analytic model, the empirical pressure- 
relation that sets their molecular fractions, their stellar and 
AGN heating mechanisms, the optical depths of their exci- 
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tation states, their assumed dependence on metallicity, and 
the normahzation factor of their CO luminosities were all 
calibrated from low-redshift measurements. These studies 
have primarily lacked a detailed model for giant molecular 
clouds (GMCs) in the ISM of high-redshift galaxies that de- 
scribes how molecular fractions, GMC sizes and densities, 
gas temperatures, and molecular level populations depend 
on halo mass, redshift, and metallicity. 

In principle , these questions can be addressed by numer- 
ical simulations. iNaravanan et al.l (|2012h has conducted sim- 
ulations of galaxies in isolation out to 2: ~ 2 with a modified 
version of the smoothed-particle hydrodynamic code GAD- 
GET2 coupled to the radiative transfer code SUNRISE to 
develop a general model for Xco, the ratio of molecular gas 
density to CO emission temperature, including physically 
reasonable prescriptions for the effect of metallicity, molecu- 
lar fraction, and gas temperature. However, their work lacks 
some of the important elements necessary for the calculation 
at 2 > 6, namely, the ability to extrapolate to the higher CO 
excitation states observable with ALMA at high-redshift, 
the consideration of an (often dominant) CMB background, 
and a cosmological context for their isolated galaxies. More- 
over, it is unclear whether the sub-grid prescriptions for ther- 
mal supernova feedback in the GADGET code are appropri- 
ate in the early universe as they are also tuned to reproduce 
observations at 2 = 0. 

Ultimately, modern simulations can still only probe a 
narrow dynamic range of size scales. Simultaneously cap- 
turing the full cosmological context of halo growth and the 
physics of photo-dissociative regions (PDRs) within GMCs 
is typically outside their grasp. Our approach is to use a 
fully analytic framework that pieces together physically- 
reasonable models at all scales. While some of our assump- 
tions will be quite idealized, they capture the necessary 
physics to achieve our goal of understanding the physics of 
CO transitions and predicting their flux in a high-redshift 
context while minimizing our reliance on low-redshift obser- 
vations. 

In building our formalism, we will work within a 
paradigm where star formatio n is a function of gas sup- 
ply and stellar feedback (e.g., Hopkins fc QuataertI I2OIII : 
ikopkins etlo] I2OIII : lOave et al.1 l2012h . We begin with a 
relation between star form ation, cold-flow accretion, and 
momentum-driven feedback (jPave et al.ll201 that is known 
to accurately describe the star formation rates and UV lumi- 
nosities of 2 > 6 LBGs, their dependence on halo mass, and 
their evolution with redshift ([Munoz 2012). We incorporate 
these inflow and outflow rates into a model for radiation- 
pressure-supported galactic disks ori ginally developed to de - 
scribe 2 ~ 2 starbursts and AGN ( Thompson et alll2005l ) 
but i mproved and exten ded to apply to higher-redshift sys- 
tems unoz fc Furlanet to 2012). This disk model provides 
radial distributions of star formation rate, gas density, and 
effective temperature for assumptions about the angular mo- 
mentum transport mechanism within the disk (^l]). We use 
these quantities to calculate the masses, sizes, and number 
of GMCs as a function of disk radius ( ij3.ip and propose a 
new method for determining molecular fractions such that 
the star formation law is obeyed at each radius given the star 
formation rates and gas densities prescribed by our frame- 
work ( i)3.2|l . On the smallest scales, we consider PDRs within 
GMCs, the geometric extent of CO-enriched gas ( i)3.3p . the 



gas temperature within clouds O, and the state of the ro- 
tational excitation levels ([JS]). In ij6l we compute the total 
CO line luminosity from the cloud and CO excitation prop- 
erties calculated in the previous sections by counting the 
number of clouds as a function of galactic radius. We com - 
pare our predictions for Xco to the Naravanan et al.] (|2012l ) 
model in !J7] and, in [JS] make predictions for the CO lumi- 
nosity function ( i]8.1|) and for foUowup observations of UV- 
selected targets ( H8.2p . In ^ we compare our results with 
those from previous stud ies, particularly from the work of 
lObreschkow et al] (|2009bl ). Finally, we summarize and con- 
clude in 



2 GALACTIC DISKS 



In iMufioz fc Furlanetto! {201^ ). we developed a one-zone 
model for the ISM of galaxies at 2 > 6. Here, we briefly 
outline the model far enough to calculate the quantities rel- 
evant for the present work, but we refer the reader to our 
previous paper for further detail s. Our formalism is based on 
work bv [Thompson et al.l (|2005l ) but improved and tailored 
speciflcally to describe z > 6 systems. 

The model assumes a rotationally-supported galactic 
disk in which marginal Toomre-instability (i.e., Q = 1) and 
vertical hydrostatic equilibrium are maintained by radiation 
pressure from stars and mechanical pressure from super- 
novae. Embedded in an isothermal halo of mass Mh at red- 
shift z whose gravity dominates the rotation, the disk has a 
gas density as a function of radius given by 



P = 



B2 



(1) 



where = y/2a/r and a is the halo velocity dispersion. 
B — 1 ~ 8.9 is the ratio of the vertical self-gravity of the 
disk to that of the halo, so that the total vertical pressure is 



Ptot ^Bph^ n^ 



(2) 



where h is the vertical scale height. The disk sound speed 
is, thus, given by Cs = B^''^ hfl. The thickness of the disk 
as a function of radius is set by the amount of gas flowing 
through it toward the center: 



h = 



M 



4:77 r p Vin 



(3) 



where vin is the inflow velocity and M is the mass inflow 
rate. Additionally, we take the maximum disk radius, R^, 
to be a fixed fraction \/y/2 of the halo virial radius with 



A = 0.05 (|Mo et al.lll99S 

The physical properties of the disk are thus determined 
by the mass and redshift of the host halo and the inflow 
of gas through the disk as a function of radius. We set the 
boundary condition on the latter so that the infall rate onto 
the edge of the disk is given by the cold flow accretion rate, 
Md, the a verage of which as a function of halo mass and 
redshift: is (jMcBride et al.ll2009l ) 



3M0/yr 



VlOio Mq 



1+2 

7 



.fb 
0.16 



(4) 

where /b is the cosmic baryon fraction. This choice pro- 
vides a cosmological context for our galaxy models fed by 
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the growth of large-scale structure. At each radius from the 
edge, star formation and outflows deplete the amount of in- 
flowing gas by an amount 2 7rr E* (1 -I- rewind) dr, where E* 
is the star formation rate surface den sity and ?7wind is the 
wind m ass-loading factor. We found in lMuiioz fc Furlanettd 
(|2012h that the pressure contributions to hydrostatic equi- 
librium of momentum transfer from ionizing photons on neu- 
tral galactic hydrogen and of supernovae are nearly equal so 
that, where pressure from stars dominates, the star forma- 
tion rate surface density is 



2ec 



(5) 



where the efficiency with which star formation converts 
baryons into radiation, e « 10""^, is only weakly depen- 
dent on IMF. Moreover, we parameterize rewind as owing to 
momentum-driven winds, 



T?wind = Vo ■ 



100 km/s 



10- 



(6) 



and set -qo = 4 to match mo del fits of the high-z UV luminos- 
ity function (iMunoj 120121') and in agreeme nt with numeri- 
cal simulations ( Oppenheimer fc DavabOOa ) . We further as- 
sume that star formation heats the gas to a temperature, T, 
given by 



o"SB T 



1 ^ 2 



(7) 



such that, in the inner regions of the disk — but still well 
outside the accretion disk of a central black hole — thermal 
gas pressure, 



p — 



pkBT 



(8) 



may contribute to the maintenance of vertical hydrostatic 
equilibrium. 

Finally, we consider two different phenomenological 
models of angular momentum transport in the disk in which 
global gravitational torques, rather than local processes, set 
the inflow velocity of gas, Vin, as a function of radius. In a 
disk with a linear spiral wave (LSW), Win is a llowed reach a 
fixed fraction of the local sound speed (e.g., iKalnaisI Il97ll : 
iLvnden-BeU fc Kalnaislll972l : lGoodmanli2003D : 



Vln =mcs, (9) 

where m is the Mach number. On the other hand, linear 
inflow due to shocks and orbit-crossings gives a velocity 
that is a constant fraction, of the local circular velocity 



tnat IS a constant traction, p, ot tne local circular velocity 
fe.g..|Papaloizou fc Pringld[l977l : iBinnev fc Tremaine|[l987l : 
iHopkins fc Quataertll201ll '): 

Win = ^r-f7. (10) 

We showed in iMunoz fc Furlanettd (|2012|) that inflow driven 
by a local a- viscosity ( Shakura fc SunvaevI 1 19731 ) cannot 
supply enough gas to the center of the disk to grow a black 
hole. Therefore, we ignore such a model in the present work. 

Combining equations IHIOI the model outlined in this 
section gives E*, pg, h, and T as functions of radius for a 
specified halo mass and redshift. Figure [1] shows the gas and 
star formation rate surface densities as functions of galactic 
disk radius in 10^'' M© halos as z = 6 for three fiducial angu- 
lar momentum transport models: LSW models with m = 0.2 
and 1 and a shocked infall model with /3 = 0.01. Note that. 
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Figure 1. Disk properties from the model described in [|2]for a 
halo mass of 10^*^ Mq at z = 6. The top panel shows the gas 
surface density as a function of radius for three different angular 
momentum transport models: LSW models with m = 0.2 (solid, 
red) and 1 (short-dashed, magenta) and nonlinear infall with /? = 
0.01 (long-dashed, blue). The horizontal line marks SSMq/pc^, 
the minimum surface density of GMCs. The lower panel plots the 
star formation surface density for the same fiducial models. 



since each of these models has the same halo mass and red- 
shift, the total gas masses and star formation rates are the 
same in each case. Results for a "high" inflow shocked model 
with P = 0.1 are similar to those for the LSW model with 
m = 1. For one of these models in which gas is transported 
to the center of the disk very quickly, the radial distribu- 
tions of galaxy properties will be more uniform throughout 
the disk. In the following sections, we will use these disk 
properties to characterize GMCs in high-redshift galaxies 
and calculate their CO line-luminosities. 



3 THE MOLECULAR GAS 

We now consider the properties of the molecular gas 
in the disks described in ij2l While we use the results 
from the previous section — which assumed strict cylindri- 
cal symmetry — we focus our attention on smaller scales 
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whe re the sas is cl u mpy a nd arranged into clouds. We adopt 
the IWolfire et alj (|2010l ) picture of multilayered photo- 
dissociative regions (PDRs) within clouds (see the schematic 
diagram in their Fig. 1) and calculate GMC properties 
(i ]3.1|) . the molecular fraction ( ii3.2|l . and the dark gas frac- 
tion ( H3.3|l as functions of r adius through the disks of high- 
redshift galaxies. Following iKrumholz et all (|2009bl ). we as- 
sume that all clouds in regions where Eg < 85 Mq pc~^ have 
constant and uniform surface densities of Ed = 85 Mq pc~^. 
Otherwise, Eci = Eg. This threshold is plotted against Eg for 
a 10^° Mq halo in the top panel of Figure [T] Note that most Substituting for Md from equation [TTI yields 
of the gas mass and star formation occur in the outer por- 
tions of the disks where the GMCs will typically be denser 
than those found locally. 



M(rci) oc r^. Thus, if all clouds have the same surface den- 
sity, then the total surface area of all clouds is proportional 
to the total mass in clouds. Moreover, the typical GMC ra- 
dius is _Rh2 = ,/h2^ -Rci- 

The distribution of clouds throughout the disk, 
dN/dlnr, is given by the ratio of the total molecular mass 
at a given galactocentric radius to the GMC mass: 



dN 2^rVH, Eg 

dlnr Mho 



2 7rr'^ Eh 



dla?^'^' ^ 7r3 G2 E2 ' 



(13) 



(14) 



3.1 Properties of Giant Molecular Clouds 

In this section, we calculate the mass, number, and visual 
extinction of GMCs in z > 6 galaxies. We assume that all 
gas in the disk fragments into atomic-molecular complexes 
with masses of order the local Jeans mass: 



Mcl = Mjcans = 



TT-* T.i 



(11) 



where B « 9.9 is defined in 3U Note that this mass includes 
both molecular and atomic gas as the fragmentation process 
does not distinguish between the two. The typical molecular 
mass within one of these clouds is given by ~ /h2 Md, 
where is the molecular fraction of the disk. This frame- 
work implicitly assumes that all gas fragments in this way 
and, in particular, that these is no smooth molecular com- 
ponent. We present our new method for determining in 
i]3.2l but for the purpose of calculating cloud properties in 
this section we give the result as 



/h2 



0.8 



[/ 0.1 Gyr 



^"cl,6 ^cl,85 '^5 



(12) 

where M^i.e = Md/lO'' Mq and Eei,85 = Eei/(85 Mq pc'^). 
The top panels of Figure [5] show the variation in cloud mass 
with fractional galactocentric radius within z — 6 galaxies 
for our fiducial angular momentum transport models over 
a wide range of halo masses. Generally, lower gas surface 
densities and increased angular velocities in the central por- 
tions of the disk cause the Jeans mass — and so the GMC 
mass — to plummet from the outer disk edge to the disk cen- 
ter. The cloud mass eventually drops bellow a solar mass so 
that stars themselves should have difficulty fragmenting out 
of the disk. While this may be unphysical and demonstrate a 
limitation of our model, hydrostatic equilibrium in these re- 
gions is mostly maintained by thermal pressure rather than 
by stellar feedback. Thus, the small cloud masses do not 
undermine out calculation of disk properties there. 

Atomic-molecular complexes have a typical radius of 
Rci ~ \/Afci/(7rEc i) such that, when Eg > 85M0pc"^, 
Rci ~ h. Following IWolfire et all (|2010l ). the average num- 
ber density within a cloud is fid = 3 Md/(47r iici A^h) = 
3 A^d/(4-Rd), where A^d is the column density of the cloud 
and hh = 2.34 x 10"^* is the proton mass. The cloud sur- 
face density, Ed, is assumed to be constant with radius, 
rd, within the cloud, and thus the average density scales as 
n oc rZ^. This means that the mass enclosed within rd is 



As shown in the central row of panels in Figure (2] the num- 
ber of GMCs increases strongly toward the center of the disk 
due to the rapid drop in cloud mass. 

The average visual extinction coefficient within the 
GMC, Av, will be important in ii3.3l when considering 
how deepl y dissociating radiation penetrates into the cloud. 
Following rWolfire et all (|20ld ). we assume, for simplicity, 
that the dust-to-gas ratio is constant within the cloud and 
that the dust opacity scales with metallicity. Normaliz- 
ing the extinction to be unity for a column density of 
1.9 X 10^^ Z'''^cm-^ gives 



Av = 



Ad Z' 



1.9 X 1021 cm- 



(15) 



where Z' — Z/Zq is the metallicity in solar units. As shown 
in Figure [2l where Eg > 85Mqpc^2^ j\^y gc Eg. Otherwise, 
Av is constant. 



3.2 The Molecular Fraction 

In this section, we consider a new approach for calculating 
the molecular gas fraction in the disk, /hj ~ M-f^^/Md, that 
attempts to recon cile two opposing views of star formation. 
On the one hand, IKrumholz et al.1 (|2009bl l express star for- 
mation as a purely local process depending critically on the 
amount of molecular hydrogen present at a given location: 

E. = Eg/H2^, (16) 
Iff 

where SFRfr/tff is the star formation efficiency per free-fall 
time divided by the free-fall time. After modeling /h2 and 
the efficiency as the fundamental physical in puts, they cal- 
culate E* as the derived quantity. In contrast. iHopkins et al.l 
(|2012h find that galactic star formation depends on a global 
competition between the infalling supply of cold gas and 
feed back and negligib ly on local molecular gas chemistry (see 
also lDave et al.|[2012l l: this latter perspective was implicitly 
assumed in the model presented in !j2l Thus, we take E* 
as the primary input. Calculating SFRff/tfi from the GMC 
properties determined in i )3.1l we solve for the required /h2 

fro m equation | 16l 

IKrumholz fc McKed (120051 ) derive SFRff /ts as 



SFRff 



-0.32 a^-- 1/4 ^3/4 



(17) 



ts 53 Myr 

where Mh2,6 = Myi^/I^^Mq, and IKrumholz et all (|2009bD 
adopt Ecore = 0.3 aud a fixed virial ratio of a^ir = 2. The 
ID Mach number of the turbulence in the GMC is given by 



© 0000 RAS, MNRAS 000, 000-000 



M,/Me=10« 



MyM„=10"' 



GMCs and C0atz>6 5 
M,/Me=10'2 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
r/Rd r/R, r/R, 



Figure 2. GMC properties for our fiducial angular momentum transport models in 10^" Mq halos at 2 = 6. The top, middle, and bottom 
panels show the GMC mass, the abundance of clouds, and the visual extinction through a cloud, respectively, as functions of galactic 
radius. As in Fig.[T] solid (red), short-dashed (magenta), and long-dashed (blue) curves correspond to LSW models with m = 0.2 and 1 
and a nonlinear infall model with /3 = 0.01, respectively. 



M = CTci/cs , where ad is the cloud velocity dispersion given 
by 



2 

0"cl = 



(18) 



5 R112 5 

Substituting for ad and = B^''^ hO,, the Mach number is 

M = 0.017 a]/;} M^/^g Ey_^5 [h/r)-\ (19) 

where aso = cr/SOkms"^. These Mach numbers are typi- 
cally supersonic for the extremely thin disks in our mod- 
els. However, they are somewhat lower throughout most of 
the disk than would expected if thermal gas pressure were 
the dominant pressure source. Substituting into equation 1 171 



with Ecorc = 0.3 and Qvir = 2, we have 
SFRff 



r-0.33 ^0.32 /. /_n0.32 



tfi 86Myr ' ^ ' 

The higher pressure in the majority of the disk cre- 
ates a somewhat higher ef Bciency than anticipated in the 
iKrumholz fc McKeel (|2005h formulation for local clouds. 

Substituting this expression into equation [16] and solv- 
ing for gives 

1.52 



/h2 



0.8 



0.33 -0.32 



[h/r 



(21) 

Equation 1211 was derived with no concern for the physical- 
ity of /h2; if the star formation rate surface density at a 
given radius is higher than can be accounted for by the 
total gas present and the efficiency specified by equation 
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Figure 3. The molecular fraction, fn^ , as a function of galactoeentric radius for our three fiducial disk models z = 6. Solid, black curves 
in the left, center, and right panels show our results from Eg. ll2l for LSW models with m = 0.2 and 1 and a nonlinear infall model with 
/3 = 0.01, respec tively. We assume host h alo masses of 10* (top), lO^O (middle), and lO^^ Mq (bottom). Dashed curves show /jj, given 
by the model of lKrumholz et"aL I ll2009d) in Eo. 1221 with c = 5. Each dashed line from top to bottom is calculated for Z' = 1, 0.3, 0.1, 
0.03, and 0.01. 



1171 then equation [21] will simply return a molecular fraction 
greater than unity. The physical implication of this scenario 
is that the clouds need to be more efficient than predicted 
by equation [20] to produce a star formation rate capable of 
supporting Q = 1 in the disk. However, except for the most 
massive galaxies, our method typically produces results for 
in a physical range. Nonetheless, to account for the pos- 
sibility, we simply set /h2 to be the maximum of equation 
1211 and unity. The result is the expression previously given 
in equation 1121 

We plot the resulting molecular fraction in Figure [3] for 
our fiducial galaxies as a function of radius and compare 
it to the values obtai ned using the radiative and chemical 
equilibrium method of lKrumholz et al] (|2009al ). In this pre- 
scription, is an analytic function only of the surface 



density of the cold gas complex, Ecomp 
clumping factor, and of metallicity: 



/h2 



1 + 



4 l + (5 



-1/5 



(22) 



where s = ln(l + 0.6 x)/(0.04 Ecomp,i ^'). X = 0.77(1 + 
3.1 z'°-^^-'), S = 0.0712 (0.1 + 0.675)"^ *, and E,omp,i = 
Scomp/(lMo pc~^)El Since the clumping factor is sensitive 



^ We have presented th e slightly updated approximation given in 
iKrumhoIz et al.l l|2009bl '). which is more accurate for small molec- 
ul ar fractions. However , we op t not to use the further refinements 
of lMcKee fc Krumhol3 l|2010r i: though applicable to lower metal- 
licities, this further method introduces additional dependences on 
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to spatial resolution, it i s difficult to determine ob servation- 
ally. However, we follow iKrumholz et alj (|2009bl ') in adopt- 
ing c = 5 for our comparison. According to equation 1221 the 
gas becomes fully molecular when Eg > 10/{cZ') Mq pc~^. 
By comparing with the top panel of Figure [TJ we can see 
that this should always be the case for galaxies hosted by 
10^° Mq halos for each of our three fiducial disk models if 
cZ' > 1. Figure |3] shows the full radial dependence of 
calculated using this method for the densities given by our 
galaxy models an d a range of meta ll icities between Z' — 0.01 
and 1. However, IKrumholz et al.| (|2009bh note that, while 
their general method is applicable to gas of arbitrary metal- 
licity, the particular approximation of equation [22] is not 
valid for Z' < 0.05. Thus, the extrapolations to lower metal- 
licities should be evaluated accordingly. 

Given its ease of use and transparent physical origins, 
equation [22] is quickly becoming a common means of cal- 
culating the molecular fraction in semi- analytic and nu- 
merical models of galaxy formation (e.g. [ Naravanan et al.l 
l2012l : IVallini et al1l2012l : iKuhlen et al.ll2012l ). However, our 
model does not provide a simple prediction for the gas- 
phase metallicity of 2; > 6 galaxies. Observational determi- 
nations are difficult because of the complexity of and degen- 
eracies within stellar population synthesis m odels in addi- 
tion to unresolved neb ular emission lines (e.g.. lDunlop et al] 
l2012l : ld"e Barros et al.|[2012i ). Numerical simulations predict 
metallicities in the range of Z' « 0.1 — 0.5, but these 
may be over-predictions especially in the f aintest galaxies 
ISalvaterra et al]|201ll : iFinlator etaLll201ll '). On the other 
hand, while our calculation of does not depend directly 
on metallicity, ou r fitting to the UV luminosity function in [J2] 
fsee lMunozll2012l) may have encoded metallicity information. 
Fixing c = 5, we ob tain rough agreement bet ween our results 
and those using the lKrumholz et al.l lj2009b[ ) method for fn^ 
by naively extrapolating down to Z' ~ 0.03, somewhat more 
met al-poor than the pre dictions of ISalvaterra et all (|201ll ) 
and iFinlator et al] (|201ll ). However, we can obtain the same 
agreement at Z' ~ 0.1 if molecular clouds in high-redshift 
galaxies are somewhat less clumpy with c ~ 2. We leave a 
more thorough exploration of the implications of these re- 
sults for the local GMC chemistry and structure to future 
work. However, in the next section, the generally low value 
of Z' will crucially limit the amount of carbon in the form 
of CO. 

In our model here, we have considered only star forma- 
tion in molecular gas. However, star formation may proceed 
in atomic gas if the metallicity is low enough that thermal 
equilibrium is reached with in a few cloud fr ee-fall times but 
chemical equilibrium is not (|Krumholj2012l ) . While this sce- 
nario typically requires metallicities below a few percent of 
solar, i.e., at the very low end of what is predicted here 
for z > 6 galaxies, such circumstances may be more eas- 
ily achieved for the decreased free-fall times of the higher 
density clouds at high-redshift. If this effect is significant in 
the galaxies we consider, our method will somewhat overes- 
timate the molecular fraction and, ultimately, the resulting 
CO luminosity. 



3.3 The Dark Fraction 

In this section, we consider how much of the carbon in 
a cloud is actually in the form of CO. Because H2 self- 
shields more effectively than CO, carbon must rely more 
heavily on dust to maintain its molecular state. This intro- 
duces a metallicity dependence into the CO fraction, fco, 
the fraction of the total cloud gas mass containing CO. 
The fraction of molecular gas mass, thus, not probed by 
CO observations — often galled the dark fraction — is /df = 
1 — .fco //h2 ■ Because the mass enclosed within a given cloud 

radius is proportional to R^,, the radius within which carbon 

1/2 

is primary in the form of CO is -Rco = . fnn ^ ci . 

In their PDR model, IWolfire et al] (|2010l ) assume that 
turbulently-generated inhomogeneities within clouds can be 
described as uniform clumps of density ric embedded within 
more difi'use smooth density gas. The clumps themselves are 
optically thin to external radiation but provide all of the 
opacity shielding the gas deeper within the cloud. The re- 
sulting fraction of the gas containing CO is given by 



In 



fco 



-4.0 



0.53 - 0.045 In f ^ ] - 0.097 In(Z') 
\ Tic cm-^ J 

(23) 

where G'q is the external radiation field in the far UV bathing 
th e GMC in unit s of the local Galactic interstellar field found 
bv lDrainel (| 19781 '). 

As discussed in IWolfire et all (|2010l ') and as is evi- 
dent from equation 1231 the dark fraction is relatively in- 
sensitive to changes in the ambient radiation field or the 
clump density and depends on metallicity primarily through 
its effect on Av- Therefore, we require only rough esti- 
mates of G'q and Uc to evaluate equation 1231 We approxi- 
mate the external radiation field as G'q ~ 0.5 eE* c^/{1.6 x 
10~^erg/s/cm2)Q We further assume Uc — (n)^^^, where 
('^)mcd ~ ''■ci \/r+"33Vf2/4 is the mass- weighted median 
density in turbulen t clouds (see discussion and citations in 
IWolfire et aIll2010^ . We note that the higher radiation field 
expected in high-redshift galaxies affects the bracketed por- 
tion of equation [23] in the opposite direction as the higher 
densities and lower metallicities. 

In Figure [4] we plot fco as a function of radius for 
our three fiducial galaxy models and show the dependence 
on cloud metallicity. Metallicity clearly has an exponential 
effect on the fraction of gas that contains CO, though the 
details depend on the ability of hydrogen gas in dense clouds 
to shield its CO, i.e., for Ay to be high despite lo w metal- 
licity. The recent work by [Naravanan et all (|2012l ) empha- 
sizes the significant influence of a lowered metallicity on 
CO emission through this mechanism, but the treatment in 



[Obreschkow et al] (|2009bh includes only a linear dependence 
on Z'. 



4 GAS TEMPERATURE 

The disk model in Sj2] gives the effective temperature of the 
star forming disk in equation [7] However, we also want to 
consider the effect on the gas of the incident CMB, whose 



the dust cross-section, the external radiation field, and the clump 
density — none of which can be usefully constrained in these galax- 
ies. 



^ Interestingly, this reproduces empirical Kennicutt relation be- 
tween star formation rate and FIR luminosity. 
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Figure 4. The fraction of all gas, /cOi in the central, CO-containing region of a cloud as a function of radius for our fiducial disk 
models in 10^'' halos at 2 = 6. As in Fig. [T] solid (red), short-dashed (magenta), and long-dashed (blue) curves correspond to LSW 
models with m = 0.2 and 1 and a nonlinear infall model with /3 = 0.01, respectively. For each disk model, lines from top to bottom were 
calculated for metallicities of Z' = 1, 0.3, 0.1, 0.03, and 0.01. In cases where fewer than four lines of a given type are shown, only the 
highest metallicities appear on the plot. 
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Figure 5. The local effective temperature of the disk, T, as a function of radius for our fiducial disk models in 10^" halos at 2 = 6. 
As in Fig.[Tl solid (red), short-dashed (magenta), and long-dashed (blue) curves correspond to LSW models with m = 0.2 and 1 and a 
nonlinear infall model with /3 = 0.01, respectively. For comparison, the horizontal, dashed line denotes the CMB temperature at 2 = 6. 
For calculating CO level populations and fluxes, we set the gas temperature to be Tg^^g = T* -|- T^j^jg . 



temperature at high redshift is much higher than at 2: = 0. 
We assume that some combination of dust and line absorp- 
tion allows the CMB to impart a significant fraction of its 
energy to the gas, energy that then contributes to the kinetic 
motions of all particles in the gas. If the entire spectral en- 
ergy of the CMB is absorbed in this way, the net effect is to 
raise the gas temperature to a value: 

Tgis =T' + Tcmb, (24) 

where Tcmb = 19.11 K (1 -f z)/7 is the temperature of the 
CMB and T is given by equation [T] In principle, this addi- 
tional energy also contributes to the thermal pressure of the 
gas in equation[51 However, since we expect more CO emis- 
sion from the larger area in the outer portions of the disk 
where radiation pressure dominates, we ignore this extra 



contribution to the thermal pressure for simplicity. While 
equation [21] includes the effect of the CMB, we have ne- 
glected the contribution from cosmic rays. This choice is also 
primarily for simplicity given the uncertainty in extrapolat- 
ing their significance to z > 6 from low redshift determina- 
tions. In Figure [5l we plot the effective temperature of the 
disk as a function of fractional disk radius for galaxies in our 
model at z = 6 and compare with the CMB temperature at 
the same redshift; the combination of these two tempera- 
tures determine the gas temperature through equation 1241 
Heating from the CMB will be particularly important for 
galaxies in low-mass halos but likely not significantly affect 
the most UV-bright systems. 

In contras t to our method in equation 1241 
iKrumholz et al.l (|201ll ) provide a physical determina- 
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tion of Tgas by considering gas in thermal equilibrium 
heated by cosmic rays and the photoelectric effect of dust 
grains and cooled by atomic and molecular line emission. 
This is the method in corporated into the simulations of 
[Narayanan et al.l(|2012l ). However, their treatment does not 
suit our purposes for two reasons. First, heating from the 
CMB, which can be neglected relatively safely at 2 ~ 2 
(their definition of "high-redshift" ) , can be a dominant 
effect at 2 > 6 and is not included in their model. Secondly, 
the cloud densities and extinctio ns in our galax i es are 
typically in a regime where the iKrumholz et al.l (|201ll ) 
model does not provide consistent solutions (i.e. the hatched 
regions of their Fig. 1). 



5 CO EXCITATION 

To calculate the optical depth and luminosity of CO transi- 
tions, we must consider the state of their level populations. 
It is not valid to equate Tgas with an excitation temperature 
if the level populations are not in thermal equilibrium. In the 
rest of this work, we consider rotational transitions J — >■ J— 1 
having rest-frame frequency fco J, where J is the rotational 
quantum number of the upper state and uco = 115.3 GHz 
is the frequency of the CO(l-O) transition. In particular, we 
are interested in the J = 1 transition, the most likely to 
be thermalized, and the J — 6 transition, since its observed 
frequency ai z = 6 falls into the lowest energy ALMA band. 

Because we have information about the structure of the 
clouds within our model galaxies, we can consider the level 
populations and optical depth of CO in each cloud directly. 
We make the simplifying assumption that the cloud is uni- 
form within a radius Rco with a density, fico = fi{< Rco), 
given by 

nco = 3 iVci/(4 Rco) « 40 cm"^ E^^f^^ Ad^l ./co • (25) 

The optical depth of the J — >■ J — 1 transition at the line 
center is 



TJ = 



gj-i 4(2 7r)3/2cr, 



nj-i Rco 1 



nj g.j-1 



n.j-1 g.j 



(26) 

where g,j = 2 J -I- 1 is the degeneracy of the Jth state, 
Aj^j-i is the Einstein-A emission coefficient for the transi- 
tion, Aj = c/Jvco, and nj is the density of CO molecules 
in the Jth s tate. We ta ke values of Aj^j-i from the online 
data files of ISchoier et al. ( 20051^ updated with recent re- 
sults from lYang et al.r ( 201ol ). The density in the J — 1 state 
is given by n,7_i = [3nco/4-Rco] [w(CO)/nH] ^co,j-i, 
where n(CO)/nH ~ 2 x 10~* Z' is the CO abundance of 
the gas and Zco,.j-i is the CO partition function in the 
J — 1 state. In thermal equilibrium, nj gj-i/nj-i gj = 
exp(— ftp Ji^co/fcb Tgas) and Zco,j~i ~ exp[—hpJ{J- 
1) vco/kh Tgas]/\/l + (fcb Tgas/ftp fcoY- Because nco ^co 
is independent of fco , the optical depth depends on metal- 
licity only as rj oc Z' in thermal equilibrium. 

Given an optical depth tj for the J — >■ J — 1 transition, 
only a fraction /3j of the line photons produced escape the 
cloud on average, where (|Drain 3 I2OIII ) 



This trapping of photons decreases the critical density neces- 
sary for thermal excitation of a transition because radiative 
excitation will supplement coUisional. This effective critical 
density is given by 



kj 



(28) 



where is the coUisional de-excitatio n rate coeffic ient 

for the transition by H2 taken from .Schoicr et al.l (|2005l . see 
footnote [3|. Emission from the CMB will have a similar ef- 
fect on the level populations, particularly in the CO(l-O) 
transition and where the line is sub-thermally populated and 
optically thin (see Appendix I A3|) . 

In practice, fico and tj span a wide range of values 
depending on metallicity, halo mass, and galactocentric ra- 
dius. Thus, certain lines will be thermal in some regimes but 
not in others. Moreover, the determinations of level popu- 
lations and optical depth are coupled due to photon trap- 
ping in thick regions. Rather than commit to particular as- 
sumptions for optical depth and t hermalization, we adopt 
the flexible approach provided by IKrumholz fc Thompson! 
(200^, which solves the coupled equations of detailed bal- 
ance assuming the CO is in chemical (even if not necessar- 
ily thermal) equilibrium. We use a modified version of their 
publicly available codcQto calculate the populations, optical 
depths, and emissivities of each line as a function of cloud 
properties (see Appendix IX)) . Of course, this code has its 
own built-in assumptions and simplifications, such as using 
a constant density throughout the cloud for the calculation 
of level populations, but it currently remains the best avail- 
able approach for addressing this problem. 

To demonstrate the variation in level populations and 
optical depth, we plot in Figure [6] the frequency-integrated 
emissivity of clouds as a function of fractional galacto- 
centric radius for a wide range of halo masses and as- 
suming our fiducial LSW model with m = 1 as an illus- 
tr ative case. We contrast results from the modified code 
of IKrumholz fc Thompsoiil (|2007l ) with what would be ob- 
tained by assuming thermal equilibrium with tj ^ 1 and 
Tj ^ 1. We further compare our results to a model with an 
extra level of complexity such that line emission is thermal- 
ized and optically thick in regions of the cloud that both 
contain CO and where nco > ncrit. For the purposes of 
this comparison, we adopt ncrit = 10'^ cm~^ for J = 1 and 
Merit = 2 X 10^ cm"'' for J = 6, i.e., the approximate val- 
ues for /?,/ = 1. We find that, in the inner regions of large 
galaxies or in smaller galaxies where the CO-containing re- 
gion is small and dense, the emission is typically thermal 
but may be both optically thin or thick. On the other hand, 
in the outskirts of the largest galaxies, the size of the CO- 
containing regions of clouds is often large enough, and the 
average density consequently low enough, for the emission 
to be subthermal, particularly for the J — 6 transition. 
However, the sub-thermally populated gas still contributes 
a significant amount to the emissivity; that is, there is more 
emission than there would be if we simply disregarded the 
portions of the cloud below the critical density. 



^ |http: / /home.strw.leidenumv.nl/~moIdata/CO.htmI| 
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Figure 6. The frequency-integrated emissivity (luminosity per unit volume of the CO-containing region) of J = 1 (top panels) and 
J = 6 (bottom panels) transition of CO as a function of fractional galactocentric radius. We set z = d and assume an LSW model with 
m = 1. The galaxies are ho sted in 10^ (left panels), 10^*^ ( middle panels), and 10^^ Mq halos. Solid (magenta) line shows results from 
our modified version of the iKrumholz fc Thompson! ll2007l ) code, while short- and long-dashed curves assume optically thin and thick 
emission, respectively, and a thermal population of states. The dotted curve shows results for thermal, optically thick emission in regions 
of the cloud where the average density exceeds ^crit — 

103 (J = 1) or 105-3 cm-3 (J = 6). 



6 CO LUMINOSITY 

6.1 Emission from Molecular Clouds 



The IKrumholz fc Thompsoiil (|200if ) code gives the 
frequency-integrated luminosity density, L'j, for each 
line. The ne/ocify-integrated luminosity of the cloud is given 
by {4:-k/3) Rqq c L'j /{J uco)- While clouds at the outskirts 
of the disk have lower emissivities than those near the 
center of the disk, they also have much larger regions that 
contain CO. By combining this cloud luminosity with the 
number of clouds as a function of radius in the galaxies (Eq. 
I13p . our model predicts a radial profile of CO flux in which 
the majority of emission is not concentrated at the center 
but rather produced in the outer portions of the disk. 

We can obtain the total luminosity escaping GMCs from 
each radius of the galaxy by simply counting the number of 
clouds in each disk radius: 

iivRlocL'j dN ^ 
3 J i^co dlnr ' 
Integrating equation [29] over radius will give the total lumi- 
nosity, Lj, for the Jth transition emitted by clouds in the 
galaxy. 

If the transitions of interest were in thermal equilibrium 
and optically thick, their line luminosities would be straight- 
forward to calculate. In this case, the velocity-integrated lu- 
minosity emitted from the surface at Rco by each cloud 
is 8 Rco Uci B(J uco,Tgs,B), where we have assumed that 



dLj 
dlnr 



S(J t'co, Tgas), the Planck function evaluated at frequency 
J vco and with brightness temperature Tgas, is constant 
over a top-hat line proflle with a turbulent width given by 
2 ad- In the Rayleigh-Jeans limit, B{J vco,T^mi) oc Tgas, 
and, if Tgas is dominated by Tcmb, then equation [22] de- 
pends on disk, CMC, and line properties as dLj/dlnr oc 
JV-Vco ^g^^^d^^"^'^- For > 85M0/pc^ this re- 
duces to dLj/d\nr oc f^o The remaining den- 
sity dependence in fco means that Lj is not proportional to 
Eg as is often assumed under these assumptions. The depen- 
dence of Lj on fco , moreover, makes it a strong function of 
metallicity. 

We use re l ations found in Appendix A of 
lObreschkow et all (|2009bl ) to convert the CO luminosities 
we calculate into quantities convenient for observers. The 
observed, velocity-integrated flux from each galaxy is 



Fco 



{l + z)Lj 



(30) 



where _Dl is the luminosity distance of the galaxy at redshift 
z. If the observing beam size is larger than the galaxy on 
the sky, then the flux in equation [30] can be expressed as a 
temperature as 



Tco 



■Fco, 



(31) 



2 /Cb ^CO ^gal 

where figai is the solid angle of the galaxy on the sky, and 
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Figure 7. The CO(l— 0) (top panels) and CO(6— 5) (bottom panels) flux as a function of halo mass for our fiducial disk models at ^ = 6. 
The top axis shows the average molecular surface density for a face-on galaxy corresponding to each halo mass. Long-dashed, solid, and 
short-dashed curves assume Z' = 0.3, 0.1, and 0.03, respectively. The noise levels in ten hour observations with JVLA and ALMA are 
denoted by the dotted lines. As in Fig. [S] the left, center, and right panels show results for LSW models with m = 0.2 and 1 and a 
nonlinear infall model with /3 = 0.01, respectively. 



f^gai = T^Rd/^A for a face-on galaxy, where Da is the an- 
gular diameter distance to the galaxy at redshift z. 

We can compare the resulting emission to the obser- 
vational thresholds of ALMA and the JVLA. Observing 
at 30 GHz in a 200kms~^ channel, the JVLA reaches a 
noise level of 42 /xjy after ten hours of observation. If we 
assume that the line profile emitted by the galaxy has a 
width of 2 cr dominated by the halo velocity dispersion, 
then the noise level on the velocity-integrated CO flux 
is 11.5 mjy kms^^ [J(75o/(l + z)]^^^. In comparison, the 
ALMA sensitivity per 200 km s^^ channel at 90 GHz in a 
typical ten hour observation is 47 ^iJy, leading to a velocity- 
integrated noise level of 7.5mJykms~^ [Jo"5o/(1 -|- z)]^^'^. 

Figure [7] shows the average J = 1 and J = 6 flux emit- 
ted by each of our model galaxies as a function of its host 
halo mass at 2 = 6. We plot results for Z' = 0.3, 0.1, and 
0.03 to cover the range of expected metallicities for these 
systems. Each of our different fiducial angular momentum 
transport models result in roughly similar CO luminosities, 
which indicates the robustness of our results. Comparing 
our results to the dotted curves denoting JVLA and ALMA 
noise levels after ten hours of observation, we do not expect 
either telescope to observe z — 6 LBGs in CO. Since it will 
not make a difference in terms of observations with current 
telescopes, we plot here only the average flux of galaxies 
hosted by the given halo mass. 

The behavior with increasing redshift is complicated by 



a competition between several factors. At fixed halo mass, 
galaxies become smaller and denser and their accretion rates 
from cold-fiows increase. However, these systems are also 
rarer, more distant, and potentially more metal-poor. More- 
over, the increasing temperature of the CMB provides addi- 
tional heating, greater effect on the level populations, and 
a higher observational background ( i]6.2|) . We leave a more 
detailed investigation of the subtleties of this redshift evo- 
lution to future work but note that the overall effect of in- 
creasing the redshift from 2 = 6 to z = 7 or 8 is minimal, 
no more than a factor of a few in either direction, and does 
not change the qualitative conclusion that the CO emission 
will be beyond the reach of modern facilities. 

The results in Figure [7] clearly show a nonlinear depen- 
dence of Fco on Z' as expected from the exponential relation 
in equation 1231 However, in the shocked infall model, most 
of the emission comes from the extreme outer portions of the 
disk where the effect of a lowered metallicity is less extreme 
(see Fig. |4]). As a result, Fco is somewhat less sensitive to 
metallicity in this model. 



6.2 CMB Subtraction 

lObreschkow et all (|2009bl ) emphasized that the CMB at 
high redshift should be considered, not just as a heating 
term, but as an observational background. The total flux 
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Figure 8. The CO(l-O) (top panel) and CO(6-5) (bottom panel) 
flux as a function of host halo mass at 2: = 6. Solid lines show 
the integrated flux produced by clouds in the galaxy, while short- 
dashed lines denote the required CMB flux to be subtracted (Eq. 
Vi'il assuming Z' = 0.1 and an LSW disk mo del with m = 0.2. 
The clo ud flux calculated using the method of lObreschkow et al.l 
ll2009bl . Eq. |35ll is plotted in the long-dashed curves (see dis- 
cussion in ijSJ. Dotted lines denote the noise levels in ten hour 
observations with JVLA and ALMA. 



at the frequency of each CO line from a galaxy with solid 
angle figai on the sky has components from both the GMCs 
within the galaxy and the CMB, which passes freely between 
clouds but may be absorbed if transmitted through clouds. 
Let S7ci,i, Ici.i, and Xi denote the solid angle, intensity, and 
fractional CMB transmission of the ith cloud numbered 1 to 
A'^. Assuming non-overlapping clouds the total flux observed 
from the galaxy is 



{Ic\,i + Xi /cmb) fici,i + Ic 



MB I f^gal ~ / Sicl.i 



while the CMB background from an equivalently sized area 
on the sky is 



Assuming a fraction e~ of CMB passing through cloud 
i is absorbed, where rj,; is the optical depth of the Jth 
transition for that cloud, then the background subtracted 
flux is 

JV 

-FCO-CMB = E/ ^ ("'^ ~ ^ /cmb] f2cl,i 

i = l 

= Fco — /cmb , (32) 

where 

/^CMB =8 7r'B(J^co,rcMB) y/?co^ci [l-e-'^'^^] 

(33)' 

The amount of CMB subtraction thus depends not just on 
the CMB itself but on the size and transparency of the 
clouds. A positive value of Fco-cmb corresponds to CO 
seen in emission, while a negative value indicates absorp- 
tion. 

Our methodology determines Fco and /cmb individu- 
ally to sufficient accuracy for our purposes but often predicts 
that the difference between the two should be much smaller 
than each of their absolute values. Consequently, /^co-cmb 
is very sensitive to the simplifying assumptions and uncer- 
tain details of our calculation, even up to the sign of the re- 
sult. The "raw" fluxes plotted earlier are therefore a clearer 
prediction of the model than subtracted versions, but we em- 
phasize that re al observations require CM B subtraction, as 
emphasized bv lObreschkow et al.l (|2009bl ). We compare our 
results for Fco and /cmb in Figure[8]assuming Z' = 0.1 and 
a fiducial LSW angular momentum transport model with 
m = 0.2 and show results using the CMB-subtracted signal 



7 THE X-FACTOR 

The CO luminosity is often evaluated through its predicted 
relationship to gas density through the quantity Xco, the 
ratio of molecular gas surface density to the J = 1 CO 
emission tempera t ure, c ommonly referred to as the X-factor. 
[Narayanan et all (|2012l ) used numerical simulations to pre- 
dict a general expression for Xco in terms of the average 
surface density of molecular gas and metallicity given by 



Xco ^ 1.3 X lO-" 



cm 



K kms 



M© pc-: 



(34) 



dlnr. 



where (Shj) is the mass-weighted surface density of molec- 
ular gas averaged over the galaxy. In the left panel of Figure 
[9l we compare this prediction to the intrinsic emission from 
all clouds in our model galaxies, taking Z' — 0.3 as a fiducial 
case. Note, however, that we define (Sna) ~ MH2,tot/A for 
our results, where A/Hj.tot is the total molecular gas mass 
and A is the face-on area of the galaxy. Each point along 
a plotted line showing our results corresponds to a differ- 
ent halo mass from 10^ Mq to 10^^ M© at 2 = 6. While the 
iNaravanan et al.l (|2012l ) model is only fit to simulated galax- 
ies for {Eh2 ) ^ 100 M0 /pc^ , we extrapolate their results via 
equation [34] below the regime in which they tested it. The 
agreement between our results and equation[34]is quite good 
given our very different methodology. 

In the right panel of Figure[9l we repeat the comparison 
with our model recalculated to exclude the CMB as a source 
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Figure 9. The X-factor as a function of average molecular fraction for face-on galaxies in our fiducial disk models at 2 = 6 with a 
metallicity of Z' = 0.3. In the left panel, we show results that include the CMB as a source of both heating and (de-)excitation. Thin 
lines correspond to results using our methodology with solid (red), short-dashed (magenta), and long-dashed (blue) curves denoting LSW 
models with m = 0.2 an d 1 and a nonlinear infall model with /3 = 0.01, respectively, as in Fig. [T] Thick curves were calculated using the 
[Narayanan et ah 1 1 I2OI2I) model (Eq. [34J but using our molecular surface densities; dashed portions indicate where we have extrapolated 
the fitting formula to surfaces densities below those actually simulated. We also note that our definition of molecular density — calculated 
as Mjj2,tot/^i where Mn^.tot is the to t al mo lecular gas mass and A is the face-on area of the galaxy — differs from the mass- weighted 
molecular densities of lNaravanan et al ] 1I2OI2I ). In the left panel, we show our results including the CMB as a source in both heating and 
(de-)excitation (but not as an observation al background), while in the right panel, we have repeated the calculation setting Tcmb = 
to compensate for the lack of CMB in the iNaravanan et aT ] l|2012r i model. 



of bot h heating; and (de -) excit ation since it is not included 
in the [Narayanan et all (I2OI2I ) model. Correcting for this 
difference, the collective agreement of our results improves 
for the largest galaxies that we consider. While there is still 
some disagreement in the smallest galaxies, we note both 
that equation [33] is likely not valid at these low densities 
and that galaxies in ha los as small as these m ay not be able 
to form anyway (e.g., iMuiioz fc Loebl I2OIII ). The demon- 
strated agreement for the largest galaxies — producers of the 
brightest CO emission, the most likely to exist, and the ones 
for which equation [34] is less of an extrapolation — validates 
our approach. 



8 RESULTS 

In ij6]we showed that the CO signal from 2: > 6 LBGs will be 
very difficult to detect even for ALMA. Moreover, Fco-cmb 
is sensitive to the details of the calculation because we ex- 
pect -Fco and Fcmb to be close to each other. Nevertheless, 
in this section, we present CMB-subtracted results from two 
different strategies for observing CO lines at z > 6: blind 
searches and foUowup of UV-selected galaxies. We calculate 
both the CO luminosity function and expected CO flux as a 
function of UV magnitude. This exercise demonstrates the 
functionality of our model, and the methods can be used in 
the future to predict observations for other molecular lines. 



8.1 Surveys for Galaxies in CO 

The ALMA field-of-view for observing CO transitions is ap- 
proximately 64 (6/ J) [(1 + z)/7] arcsec on a side, which can 
cover the 4.6 arcmin'^ WFC3 IR field-of-view in about four 
pointings. It will be interesting to compare the luminosity 
function of molecular line emission at high-redshift to those 
compiled from rest-frame UV data. 

After calculating the average CMB-subtracted CO lu- 
minosity, I/qo-cmb, as a function of halo mass, we use the 
[Sheth fc Tormen[ (|2002l ) halo mass function to produce the 
CO luminosity function. In principle, this calculation should 
include the scatter in luminosity for a given halo mass that 
results from variations in the galactic accretion rate from 
galaxy-to-galaxy. These fiuctuations serve t o flatten the lu- 
minosity function j ust as it does in the UV (jMufioz fc Loebl 
[201ll : [Mufio3[2012l ). However, since we do not expect these 
galaxies to be detectable in CO with current instruments 
with or without this scatter, we simplify the treatment and 
calculate the luminosity function by assuming a single value 
of the luminosity for each halo mass. Also for simplicity, we 
further ignore any cut-off of galaxy formation in low-mass 
halos. 

We show results for the CO luminosity function at 2 = 6 
in Figure [10] assuming Z' — 0.1 and a LSW angular mo- 
mentum transport model with m — 0.2, but, as shown in 
Figure [7] different disk models would yield similar results. 
We present the calculation for both the J = 1 and J = 6 
transitions. In the particular fiducial case plotted here, both 
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log,o Lco-cMB/(''y km/s Mpc^) 

Figure 10. The CO luminosity function a,t z = 6. Solid and 
dashed curves show results for the J = 1 and J = 6 tran- 
sitions, respectively, assuming Z' = 0.1 and a fiducial LSW 
d isk with m = 0.2. Our p redictions are well-below both those 
o fObreschkow et al.l ll2009bt) and the expected observational sen- 
sitivities. 

lines are seen in emission over the depicted luminosity range. 
As expected from Figures [7] and |8] our l uminosity function 
predict ions are well below both those of lObreschkow et al.l 
l|2009br i and the expected observational sensitivities. 

8.2 Followup of Existing LBG Candidates 

While telescopes like ALMA will not be as efficient as WFC3 
in discovering new, large samples of galaxy candidates at 
z > 6 based on their CO line emission, they may be use- 
ful for the followup of UV-selected targets. In this section, 
we calculate the expected CO flux as a function of Afuv, 
the absolute AB magnitude in the rest-frame UV. As with 
the luminosity function, the scatter in both CO luminosity 
and Muv at fixed host halo mass enters into the calculation 
jMunoz fc Loebl l201ll : iMunoj |2012| : [ Munoz &i Furlanetto! 
|2012| ). These fluctuations are correlated and lead to less scat- 
ter in CO flux at fixed Afuv than at fixed host halo mass. 
However, given the difficulty in observing CO, we again per- 
form a simplified computation in which every halo hosts a 
galaxy of the same average Muv and -Fco-cmb. 

Figure [11] shows the average CO flux expected from 
z = 6 galaxies as a function of absolute magnitude in the 
rest-frame UV assuming Z' — 0.1 and a LSW angular mo- 
mentum transport model with m — 0.2. Again, as in Figure 
[T] different disk models yield similar results. We have com- 
pared these results to the noise levels for measurements of 
the J = 1 transition with JVLA and to the approximate 
curre nt detection limits in the UDF12 UV lum inosity func- 
tion (jSchenker et al.ll2012l : iMcLure et al] 120121 '). Given the 
fainter flux in the J = 6 transition and the increased noise 
level observing at higher frequencies, the prospects for de- 
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Figure 11. CO flux as a function of absolute UV magnitude at 
z = 6. As in Fig. 1101 solid and dashed curves show results for 
the J = 1 and J = 6 transitions, respectively, assuming Z' = O.l 
and a fiducial LSW disk model with m = 0.2. Dotted lines indi- 
cate noise levels from JVLA observations of the J = 1 transition 
integrated for 10, 100, and 1000 hours. Shaded regions indicate 
galaxies that would cither be too faint in the UV for detection in 
UDF12 or too rare to be found in a UDF-sized area on the sky. 



tecting this line are even worse. Clearly, galaxies are more 
easily found in the UV than in CO. Only the most UV- 
bright galaxies will have detectable CO counterparts and 
only after ~ 1000 hours of integration. Stacking the 10-hour 
signal from ~ 100 such sources could achieve an equivalent 
detection if the redshift of each object is precisely known. 
Yet, these objects are also the rarest on the sky — no objects 
at 2 = 6 br ighter than M^y ~ —22 in the UV are found 
m the UDF (H ouwens et al.ll2006l ) — providing an additional 
obstacle to the stacking of large samples. 



9 COMPARISON WITH PREVIOUS WORK 

In this section, we discuss the rela tionship betwe en our work 
and recent efforts by lobr cschko w et al.l (|2009al ). These au- 
thors modeled the CO luminosity by combining a semi- 
analytic model for galax y formation (|Croton et al.l l2006l : 
iDe Lucia fc Blaizot|[2007l 'l with simple prescriptions extrap- 
olated from local empirical relations. Our formalism repro- 
duces the normalization constant for their CO luminosity 
under their set of assumptions but predicts a much lower 
value under the physical conditions expected in the high- 

reds hift ISM. 

lObreschkow et al.l (|2009bt ) give the CMB-subtracted lu- 
minosity as a function of molecular mass, metallicity, and 
excitation temperature as 

= Mn^.tot ^It^ [l-exp(-B'rj)] l{J,T,^,,z), (35) 
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where 



TcMB gives 



exp 



(^)-l -P(^) 



(36) 

Tj = 7.2 Tc exp ( — ^ , ^ ) smh ( ^ , ^ — ) , (37) 



2 /ck 



2 /i^b T-^c. 



and fc = 8 X 10"* Wkg"^ and Tc = 2 are determined us- 
ing fits to the CO luminosity function with Z' « 0.01 and 
spectral energy distributions of local galaxies. B' and e de- 
scribe the degree of cloud overlap and the amount of smooth 
molecular gas not in clouds, respectively. Since they find that 
these latter two eflects are not important at low redshift and 
roughly cancel out at high redshift, we set B' = e = 1 here 
and note that this choice also simplifies the comparison to 
our results since we ignore bot h eflects in our model. 

lObreschkow et all (|2009bl '). thus, assumes that the CO 
line emission is well-approximated by the emission from 
thermalized gas at a slightly under-estimated optical depth 
or temperature and expresses the excitation temperature of 
the CO in their model galaxies as 



CMB 



-To 



AGN) 



(38) 



where To, Tsb, and Tagn are the contributions from "reg- 
ular" star formation, star bursts, and AGN, respectively. 
Since their semi-analytic models don't provide information 
about the radial distribution of a galaxy's properties, all 
quantities are averages over whole galaxies. Since the star- 
burst and AGN components are low for typical ^ > 6 galax- 
ies, the excitation temperature is given almost exactly by 
considering only the first two terms in eQuation l38l and, from 
calibrations to observations of local systems, they adopt 
To = 17 K. 

The merging of relevant temperatures in equation [38] 
is similar to the way we combine the temperature of the 
disk with that of the CMB in equation [T] However, the 
contribution from "regular" star formation in our model is 
far from constant as is To- In the most massive galaxies 
of our model, the contribution from non-bursting star for- 
mation provides the dominant heating mechanism for CO. 
A more important diflerence with our method, however, is 
that we clearly find that CO lines are very sub-thermally 
populated in many regimes, particularly the J = 6 transi- 
tion (see FigO and employ th e escape probability code of 
iKrumholz fc ThompsonI (|2007l ') to account for this. This re- 
sults in significantly less CO flux that would otherwise be 
predicted. 

Since CMB-subtraction is not an issue for the local 
galaxies from which their model parameters are calibrated, 
we can obtain the luminosity produced by molecular gas in 
the lObreschkow et al.l (|2009bl ) model by ignoring the second 
term of eauation l36l We can further simplify equation 1351 bv 
noting that To <^ Tomb at z = 6 and that the resulting 
Tj S> 1 for both J — 1 and J = 6. This yields 



k Z' Afn, 



J" 



exp 



(39) 



By contrast, using our method to integrate the luminosity 
from all clouds as a function of galactocentric radius under 
the assumption of thermal, optically thick emission at Texc = 



Lj = 



16 TT hp Vqq I fco o-ci \ M] 



.fua 



JH2,tot 



Z' expf^i^Ul' 

yl'b J CMB / 

(40) 

where {fco ''"ci//h2 Ed) is the molecular-mass-weighted av- 
erage of fco o"ci/.fH2 Eel throughout the galaxy. Equating 
equations 1391 and 1401 gives a value for k of 



k = 6.15x10"* Wkg"^ 



0.01 
Z' 



/co//h2 85Mqpc- 



0"c 



0.5 



1 km 
(41) 

Thus, ou r methodology can r eproduc e the normalization fac- 
tor that lObreschkow et ahl (|2009tJ ) determine empirically 
under the same set of assumptions that they adopt. How- 
ever, our model also provides results when those assump- 
tions are relaxed and, in particular, allows for a realistic ex- 
trapolation to the physical conditions at high redshift. For 
metallicities in the predicted range of Z' = 0.1-0.3, k can be 
orders of magnitude lower, depe nding on halo mass and dis k 
model, than the value used bv lObreschkow et al.l (|2009bl ). 
Additionally, the higher cloud surface densities at z > 6 im- 
ply smaller, more compact clouds than observed locally and, 
consequently, less emitting surface area for the same cloud 
mass. 

A fin al important differ e nce be tween our methodology 
at that of lObreschkow et al.l (|2009bl ) is in the calculation of 
molecular gas m asses. Our metho d , descr ibed in i)3.2l is con- 
sistent with the IKrumholz et al.l (|2009al l model, which ar- 
gues that density is the most ph ysical quantity affe c ting th e 
molecular fraction. In contrast. lObreschkow et al.l (l2009bl ). 



along with several other studies (e.g., Lagos et al.ll2012l '). rely 
on the relationship between molecular fraction and mid- 
plane pressure pressure empirically determined from local 
galaxies. This relation gives a ratio of molecular to atomic 
gas of 



Rn 



(42) 



where -Rmoi = (/h2^ ~ I) ^ , P is the mid-plane pressure, a — 
0.8, and = 2.35 x 10"^^ Pa (|Lerov et al.ll2008l ). Molecu- 
lar fractions were determined from CO( l-O) measurements 
with an assumed value of the X-factor. lObreschkow et al] 
(|2009ah further assume exponential disk density profiles in 
their model galaxies and a relationship between this den- 
sity and mid-plane pressure. While this treatment may be 
sufficient for their calculations at z = 0, understanding the 
underlying physical process responsible for the atomic to 
molecular transition is critical if the relations are extrapo- 
lated to environments very different from the ones in which 
they were calibrated. Using the mid-plane pressure in our 
disk model (Eq. [2} , we find that the molecular fraction at 
z — 6 calculated from equation [42] is essentially unity at all 
radii. That is, the high densities and thin disks of our galax- 
ies result in gas that is completely molecular for all models 
we consider. With molecular fractions this high, the star for- 
mation efficiencies we calculate in GMCs would imply star 
formation rates and UV luminosities somewhat higher than 
observed. Given the relationship between molecular mass 
and CO luminosity, our lower molecular fractions will also 
result in less CO fiux. 

Figure [8] compares results for the CO fiux produced 
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by clouds in 2 = 6 galaxies using equation [35] to ours as- 
suming Z' = 0.1 and a fiducial LSW angular momentum 
transport model with m = 0.2. In evaluating equation 1351 
we assume that the molecular and gas fractions are both 
uni ty so that Mn^.tot = fh A fh. In contrast to our results, 
the lObreschkow et al.l (|2009al ) method predicts significantly 
more CO emission, which should be detectable with current 
observatories. 



10 CONCLUSIONS 
In|M 

uiioz fc Furlanettcl (|2012l ') we improved and tailored a 
one-zone ISM model, which had already been found to de- 
scribe lower redshift systems, for use at 2; > 6. In the present 
work, we used this specialized formalism to calculate the 
properties of GMCs and predict their CO line fluxes in a 
way that realistically considers the physical conditions ex- 
pected in Lyman-break galaxies at very high redshift. 

Comparing CMC properties for three different fidu- 
cial angular momentum transport mechanisms — linear spi- 
ral waves (LSWs) with m = 0.2 and m = 1 and nonlin- 
ear infall with /? — 0.01 — we found differences in the radial 
distributions of cloud mass, number, and visual extinction 
among our models. In particular, the properties for clouds 
in the m — 1 LSW model are more constant with galac- 
tocentric radius than for other models because of a higher 
velocity of centrally-flowing gas. We then proposed a new 
method for determining /hj , taking the star formation rate 
surface density, the gas surface density, and the star for- 
mation efficiency in clouds as fundamental inputs. Compar- 
ing our results with the chemical equilibrium calculation of 
iKrumholz et all (j2009bl ) suggests either that the metallicity 
of z > 6 galaxies is of order a couple percent of solar or that 
the clouds in these systems are somewhat smoother than at 
z — 0. We then use d a simple PDR model of clouds from 
IWolfire et all (|2010l ) to investigate the effect of metallicity 
on the region within the cloud capable of maintaining its 
carbon in the form of CO. 

We set the gas temperature to be a simple combina- 
tion of the local effective temperature of the disk and the 
CMB temperature at the given redshift. Using this tempera- 
ture and taking into account the higher average densities in 
the central, CO-containi ng regions of clouds, we emplo yed 
a modified version of the IKrumholz fc ThompsonI (|2007l ) es- 
cape probability code to calculate the CO emissivity with 
no a priori assumptions on the level populations or optical 
depths of the transitions. We computed the resulting fiux 
from the entire galaxy by simply counting up the number of 
clouds at each galactic radius. 

Comparing our resu lts to the universal mo del for the 
CO X- factor predicted bv lNaravanan et al.l (|2012h . we found 
very good agreement in regions that those authors sim- 
ulated when we adjusted the CMB temperature to com- 
pensate for their different methodology. However, we pre- 
dict much less emission from clouds than the semi-analytic 
method of lObreschkow et al.l (|2009bl ). This disparity is due 
to our improved treatments of molecular gas fractions, the 
dependance of CO-dissociation on metallicity, and CO level 
populations and optical depths, but we showed that our 
model can reproduce their locally-calibrated normalization 



constant for CO emission if we adopt their set of assump- 
tions. 

Our final determination of the CMB-subtracted, 
velocity-integrated CO fiux as a function of halo mass and 
redshift depends somewhat on metallicity but is fairly ro- 
bust to the choice of angular momentum transport model 
in the galactic disk. We predict that the CO signal will be 
very difficult to observe either with existing or currently- 
planned fa£ilitieS;_^n_aB_add^^ contrast to the conclu- 
sions of lObreschkow et all (|2009bl '). we also find that z > 6 
galaxies are likely absorbers of CO rather than emitters. 

Blind searches for new populations of high-redshift 
galaxies will be particularly inefficient given the large 
amount of integration time required to detect even the 
rarest objects. However, targeted studies may prove fruit- 
ful if a significant amount of observing time on the JVLA 
(~ 1000 hours) is dedicated to extremely UV-bright candi- 
dates. Because of this, shallow, wide-area observat ions in the 
rest-f rame UV, suc h as those by the HIPPIES (|Yan et al.l 
I2OIOI ) and BoRG jTrenti et all I2OI j ) prograr ns, or high- 
magn ification lensing surveys like CLASH (jZitrin et al.l 
I2OIII ) will be the most useful sources of candidates suitable 
for foUowup. 

Unfortunately, our results also imply that the integrated 
CO luminosity of the entire high-redshift galaxy population 
is much smaller than a ssumed by pr eviou s efforts. In particu - 
lar. iLidz et al] (|201ll '). ICariiiil (| 20111 ) . and I Gong et"al] (|201ll ) 
anticipated that "intensity mapping" of CO lines from unre- 
solved galaxies would provide a powerful probe of the physics 
of star formation and the intergalactic medium during that 
era. However, they all used local calibrations of the CO lu- 
minosity per unit star formatio n. Here we have shown tha t 
such calibrations, as also used bv lObreschkow et all (|2009bl ). 
greatly overpredict the luminosities. Our results suggest that 
these earlier predictions are orders of magnitude too opti- 
mistic, making intensity mapping in CO lines extraordinar- 
ily difficult. Fortunately, ALMA will be able to test our pre- 
dictions in the near future. 

In forthcoming work, we will continue to probe the ob- 
servational signatures of this model and capitalize on its 
ability to calculate CMC properties and molecular line emis- 
sion. In particular, given the dissociated state of much of the 
carbon in low metallicity GMCs, we will apply our formal- 
ism to the problem of CII emission lines from high-redshift 
galaxies. These lines have already be en observed out to z ~ 5 
in bright submillimeter galaxies (e.g. lWagg et al ] |2012l V Pre- 
dictions for more typical galaxies out to even higher redshifts 
will allow us to explore the earliest era of galaxy formation 
in the universe. 
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APPENDIX A: MODIFICATIONS TO THE 
KRUMHOLZ & THOMPSON CODE 

In this section, we document the clianges we made 
to the publicly- available e scape-probability code of 
iKrumholz fc Thompsoiil (|2007l '). 



Al Cloud Velocity Dispersion 

The original version of the code assumes that thermal gas 
pressure in the cloud maintains equilibrium against gravity 
so that the sound speed is specified by the gas temperature. 
Consequently, the input gas temperature and cloud Mach 
number determine the cloud velocity dispersion. However, 
in our disk model, the sound speed is determined primarily 
through radiation pressure and mechanical pressure from 
supernovae rather than from gas pressure (see !j2]). Thus, 
we calculate the cloud velocity dispersion and Mach num- 
ber from equations [18] and 1191 respectively, removing the 
dependence on the gas temperature. 



A2 J=l-0 Optical Depth 

The original version of the code takes the optical depth of 
the J = 1 transition as an input and determines rj for each 
molecular line as a fraction of rj=i. As such, it only does 
part of the calculation for tj before using rj=i as a nor- 
malization constant. However, the J — 1 transition is not 
necessarily in thermal equilibrium making its optical depth 
difficult to calculate. We alter the code to compute tj with- 
out the normalization constant, which requires introducing 
the radius of the CO-containing region of the cloud as an 
additional input. 



A3 CMB 

The original version of the code ignores the ability of the 
CMB to radiatively excite or de-excite levels of CO in a 
similar way to radiation from the CO itself. The absorption 
rate of CMB photons per particle is 

TcMB = Pj Bj-i^j B{J vco,Tcmb) 

^ n _9J Aj,J-l 

5J-1 e''p-'''co/fcBTcMB _ 1' ^ ' 

where Bj-i^j is the Einstein-B coefliicient for absorption 
and hp and fee are the Planck and Boltzmann constants, re- 
spectively. In an optically thin region (i.e., 13 j — 1), the 
CMB absorption rate is more significant for lower-order 
transitions. However, because an increase in a particular 
level necessitates a decrease elsewhere in the distribution, 
all levels will be somewhat affected. To modify the code, 
we note that the excitation rate per particle due to emis- 
sion from CO is (3j Aj^j-i, where the factor of /3j accounts 
for line photons that escape the cloud without contributing 
to the excitation. Therefore, we simply multiply Pj Aj^j^i, 
where it appears in the code, by an additional factor of 
(l + [s.//5J-i]/[e'''"^''°°'"'^^'="^-l]). Ultimately, this mod- 
ification reduces the luminosity in the lines only by, at most, 
a factor of a few. 
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